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Abstract 

The e~-e'^ bound state spectrum of QED3 is investigated in the quenched ladder 
approximation to the homogeneous Bethe-Salpeter equation with fermion propagators 
from a Rainbow approximation Schwinger-Dyson solution with variable fermion mass. 
A detailed analysis of the analytic structure of the fermion propagator is performed 
so as to test the appropriateness of the methods employed. The large fermion mass 
limit of the Bethe-Salpeter equation is also considered, including a derivation of the 
Schrodinger equation, and comparisons made with existing non-relativistic calcula- 
tions. 

PACS NUMBERS: ll.lO.S, ll.lO.K 
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1 Introduction 



The similarities between Quantum Electrodynamics in three space-time dimensions 
(QED3) and Quantum Chromodynamics in four space-time dimensions (QCD4) and 
the simplicity of the theory make QED3 attractive for the study of non-perturbative 
methods. QED3 is an abelian theory and provides a logarithmic confining e~-e'^ 
potential 

Our approach to positronium states in QED3 is via a solution to the homogenous 
Bethe-Salpeter equation with fermion propagator input from the Schwinger-Dyson 
equation. The full Schwinger-Dyson and Bethe-Salpeter equations are intractable. 
Here we consider a solvable system of integral equations within the quenched, ladder 
approximation. This crude truncation of the full equations does break gauge covariance 
but has very attractive features and has been employed extensively in QCD4 spectrum 
calculations |^ . 

This study continues on from a previous study Q which uses a four-component 
fermion version of QED3. In this version, the massless case exhibits a chiral-like U{2) 
symmetry broken into a U{1) x U(l) symmetry by the generation of a dynamical 
fermion mass, resulting in a doublet of Goldstone bosons. This pion-like solution is 
important for drawing similarities between QED3 and QCD4. The four component 
version of QED3 is also preferred to the two component version because the Dirac 
action in the two component version is not parity invariant for massive fermions. 
QCD4 is parity invariant and we aim to have as much in common with that theory as 
possible. 

The previous work was restricted to zero bare fermion mass, while in this study 
the bare mass is increased from zero to large values in order to compare with results in 
the non-relativistic limit. This study also takes a closer look at the choice of fermion 
propagator input. Knowledge of the analytic properties of the fermion propagator 
is important for determining the approximation's ability to provide confinement and 
whether or not any singularities will interfere with a Bethe-Salpeter solution. Based on 
the work of Maris (H, ^ the occurrence of mass-like complex singularities is expected 
which have the potential to influence our calculations. 

In section 2 we look at the Bethe-Salpeter and Schwinger-Dyson approximations 
used in this work and the method used to find the bound state masses. A brief review 
of transformation properties in QED3 is given in the appendix. These transformation 
properties are of vital importance for an understanding of the structure of the e~-e^ 
vertex function and the classification of the bound states. Section 3 describes the 
non-relativistic limit and the connection between the Bethe-Salpeter and Schrodinger 
equations for QED3. 

In section 4 the approximation to the fermion propagator is detailed. The structure 
of the propagators will be analysed in the complex plane where we attempt to locate the 
expected mass-like singularities. In section 5 the Bethe-Salpeter solutions are reported 
and comparisons are made with non-relativistic limit calculations. The results are 
discussed and conclusions given in section 6. 

2 Solving the Bethe-Salpeter Equation 

The Bethe-Salpeter (BS) kernel for this work is a simple one-photon exchange (ladder 
approximation) which is a commonly used starting point. For convenience we use the 
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quenched approximation, work in Feynman gauge and work only with the Euchdean 
metric. Fig. 1 shows the Bethe-Salpeter equation in the quenched ladder approxima- 
tion. The corresponding integral equation is 

T{p, P) =-e^J ^ D{p - q)^,S{ \ P + q)T{q, P)S{- \ P + q)^,, (2.1) 

where T{p,P) is the one fermion irreducible positronium-fermion-antifermion vertex 
with external legs amputated. The photon propagator D{p — q) in Feynman gauge is 
l/(p — g)^. The fermion propagator S is the solution to a truncated Schwinger-Dyson 
(SD) equation. A fermion propagator has been chosen which supports spontaneous 
mass generation necessary for the formation of the Goldstone bosons. The truncated 
SD equation for a fermion of bare mass m is 

np) = S{p)-^-{i^ + m)=e^ j ^D{p-q)^,SF{q)^^. (2.2) 

This approximation is the quenched, rainbow approximation named so because the 
photon propagator has been replaced by the bare photon propagator and the vertex 
function F has been replaced by the bare vertex 7, resulting in a series of Feynman 
diagrams which resemble rainbows. In the quenched approximation the SD and BS 
equations can be recast in terms of a dimensionless momentum p/e^ and bare fermion 
mass m/e^. From here on we work in dimensionless units and set = 1. 

We use either of the two following equivalent representations of the fermion prop- 
agator, 

S{p) = -ii^av{p^)+as{p^) (2.3) 
or ^ 

^^^^ " i Mip"^) + Bip"^)' ^^'^^ 

The generation of a dynamical fermion mass and the breaking of chiral symmetry is 
signalled (in the massless limit) by non-zero B{p'^). The vector and scalar parts ay 
and as of the propagator are related to the functions A and B simply by dividing 
these functions {A, B) by a quantity p'^A?{p^) + B'^{p^). 

Note that a substitution of this fermion propagator into the Ward-Takahashi iden- 
tity shows that the bare vertex approximation breaks gauge covariance. However, this 
model is simple and does meet the requirement that the appropriate Goldstone bosons 
are formed . It is not difficult to derive a zero mass solution to our BSE analytically. 
A vertex proportional to the matrix 74 or 75, defined in the appendix, will reduce the 
quenched ladder BSE to the quenched ladder (rainbow) SDE in the case of zero bound 
state mass thus forming a doublet of massless states. According to the terminology 
used in the appendix this is an axi-scalar doublet. These solutions will be seen in 
section 5. 

Once the photon and fermion propagators are supplied, the BS equation can be 
written as a set of numerically tractable integral equations. To do this, we write 
the BS amplitude F in its most general form consistent with the parity and charge 
conjugation of the required bound state, and then project out the coefficient functions 
for the individual Dirac components. It is convenient to work in the rest frame of the 
bound state by setting P^ = (0, 0,iM). Then the scalar and axi-scalar vertices given 
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in the appendix by eqs. ( |A.7| ) and ( A.9| ) can be written as 



r^(g, P) = fiqs, |q|; M) - ^t/fe, |q| ; m) 



+ ^T^lA5V{q3, |q| ; M) + ^73^^(93, |q|; Af), (2.5) 

|q| 



r^'{l,P) = j/(g3,|q|;M)-^(^;j;j jf/(g3,|q|;M) 

+^ ( ^4 ) n93,|q|;M)-73 ^ H^(g3,|q|;M), (2.6) 

where the index j takes on values 1 and 2 only, |q| = {qf + gl)^ and q"*" = (— q'2;'?i)- 
The pseudoscalar and axi-pseudoscalar vertices are obtained from the scalar and axi- 
scalar vertices by multiplication by the matrix 745. 

It is found that the same coupled integral equations result when a vertex is multi- 
plied by the matrix 745 and so (scalar, pseudoscalar) and (axi-scalar, axi-pseudoscalar) 
form two pairs of degenerate states. 

The four equations derived from the BSE, after some manipulation including an 
angular integration, are 

2 roc roo 
(27r)^ J-00 Jo (q,2_^2)2 

{Tfffiq) + TfuU{q) + TfvV{q) + TfwW{q)) 
^iP) = TTT^ / ^93 / |q|rf|q| _ '1 X 



{Tuff{q) + TuuU{q) + TuvV{q) + TuwW{q)) 

^(P) = TTT^ / / |q|c^|q| — — —r- X 



(2^)2 Jo /3(a2-/32)l 

(Tvffiq) + TvuU{q) + TvvV{q) + TywWiq)) 

^(P) = 7^ / ^^^3 / h\d\ci\ — — T X 

{Ziry J-00 Jo (q,2_^2)2 

(Twffiq) + TwuU{q) + TwvV{q) + TwwW{q)), (2.7) 



where 



a = {p3-q3r + \P\ + |qr, /3 = -2|p||q|. 
Now define the momentum Q by 

1 



=ql + \q\^ -^M-' + iMq^, (2.8) 

and use the abbreviations ay = o"v(Q^) and as = (^s{Q'^) for use in the definition of 
the functions Tff,Tfu, . . . which are analytic functions of ^3, |q|, and M. The diagonal 
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T's are given by 

Tff = {^M' + qi + \q\^)\av\'T\^s\', 

Tuu = {jM^ + qi-\cif)Wv\^±\as\\ 

Tvv = i^M^ + qi + \q\^)\av\^±\as\^ 

Tww = -(^M2 + g|-|q|2)|ay|2±|asP, (2.9) 

where the upper sign apphes to the scalar equations and the lower sign to the axi-scalar 
equations. The ofF-diagonal T's are, for the scalar positronium states: 

■■ {o-ycrs + as(Tv)\q.\, 

: -M\q\\av\\ 

■ -[U^v^s + (^*s(^v)M + iq3{a^as-<T*sav)], 

■ 2g3|q||o-yp, 



Tfu 


= ^Uf 


TfV 


= Tvf 


Tfw 


= Twf 


Tuv 


= Tvu 


Tuw 


= 


Tvw 


= Twv 



and for the axi-scalar states: 



Tfu = 


Tuf 


TfV = 


-Tvf 


Tfw = 


-Twf 


Tuv = 


Tvu 


Tuw = 


-Twu 


Tvw = 


-Twv 



Kc^V^rs - cr*sav)\q\, 
M\q\\av\^ 

-i{a^as - <T*s(7v)q3 - U^v^s + cr*s(^v)M, 
-[^{(^Ws - cr*gav)M - qsiayus + cr*s(^v)], 
-2g3|q||(Typ, 
{al,as + (J*sav)\q\. 



(2.11) 



This is the same set of equations solved in Ref. Q with only the fermion propagator 
input altered. The bare fermion mass m only comes into the calculation through this 
input. 

The solution to the BSE involves iteration of the coupled integral equations in 
Eq. ( |2.7| ). These equations may be rewritten as 

f(|p|,P3;M)= J dqs J d|q|iC(|p|,p3;|q|,g3;M)f(|q|,(?3;M), (2.12) 

where f = (/, U, V, W)^ . For each symmetry case and each fermion mass this is solved 
as an eigenvalue problem of the form 

J dq K{p, q; M)i{q) = A(M)f (p), (2.13) 

for a given test mass M. This is repeated for different test bound state masses until 
an eigenvalue A(M) = 1 is obtained. 



3 Non-Relativistic Limit 

We consider now the non-relativistic limit m ^ cx) of our BS formalism in order 
to enable comparisons with existing numerical calculations || of the Schrodinger 
equation for QED3, and with the large m solution of Eq. ( p.TD 
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The Schrodinger equation with a confining logarithmic potential is an interesting 
problem in its own right. Initially one is faced with the problem of setting the scale 
of the potential, or equivalently, setting the zero of energy of the confined bound 
states. A solution to this problem was proposed by Sen and Cornwall in 
terms of cancellation of infrared divergences in perturbation theory. They introduce 
a regulating photon mass /x in order to set the potential as the 2-dimensional Fourier 
transform of the photon propagator l/{k'^ + A*^)? leading to a potential proportional 
to ln(/ir). They further interpret the sum of the bare fermion mass and the fermion 
self energy evaluated at the bare fermion mass shell as a renormalised fermion mass, 
leading to a mass renormalisation 6m oc ln(m///). The logarithmic divergences in the 
photon potential and the fermion self energy then conspire to cancel leaving a finite 
positronium mass. 

The first numerical treatment of the Schrodinger equation for QED3 using this 
line of argument was carried out by Yung and Hamer |^]. In a subsequent, improved 
calculation by Tam, Hamer and Yung [^, the formalism was shown to be consistent 
with an analysis of QED3 from the point of view of discrete light cone quantisation. 
Their resulting expression for the bound state energy, obtained as a solution to the 
differential equation 

(--V^ + ^{C + ln(mr))l 0(r) = {E - 2m)0(r), (3.1) 

where C is Euler's constant, is given in terms of the bare fermion mass m as 

E = 2m+ —lnm+ — ( \--ln-\ (3.2) 
Air 27r V 2 Try ^ ^ 

The lightest s-wave positronium state and first exited state are given by Aq = 1.7969 
and Ai = 2.9316 respectively. The first five states are provided in ref. Q. 

Here we present a treatment of the non-relativistic limit of the QED3 positronium 
spectrum in terms of our SD-BS equation formalism. We begin with the fermion 
propagator in the limit m ^ oo. For large fermion mass we expect the residual effect 
of the chiral symmetry breaking contribution to the fermion self energy to be small 
compared with contribution from the perturbative loop expansion. We shall therefore 
assume to begin with that the self energy is reasonably well approximated by the 
one- loop result. The validity of this approximation for spacelike momenta will be 
demonstrated numerically in the next section. 

The 1-loop fermion self energy, with the functions A and B defined in Eq. ( |2.4| ) , is 
given by 

A = l + (3.3) 



where 



and 



S^(x2 



8ttx^ 



m 



( \ — 

arccos 



2x \ 1 + 



(3.4) 



r>2 



B = m I 1 + ^^i^ I , (3.5) 
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where 



Sirx 



■ arccos 



1 + x2 



(3.6) 



This result is valid for (Euclidean) spacelike momenta > 0. An analytic continuation 
of the S functions valid for < m^, or |x| < 1, is 



1 



Sttzx 



Eb(x2) 



\l — ix J X 



3 , f I + ix 
m 



Sirix 



1 



IX 



(3.7) 



(3. 



Note that this representation exposes a logarithmic infinity in the self energy at the 
bare fermion mass pole = —m?. This is the infrared divergence in the renormalised 
fermion self energy as defined by Sen |1C] referred to above. However, in our formalism, 
this singularity does not lead to a pole in the propagator functions ay and as defined in 
Eq. (^]^) , which would signal the propagation of a free fermion [^] , but a logarithmic 
zero. 



Using Eqs. (3^) and ( |3.5| ) we obtain 



cry(p^) 



1 + ^ 



^1 + ^) +2(1 + 



2m / m 



(3.9) 



0-5 (/ 



where we have defined 



1 + 



m 



' m 



) +2(1 + 



o 9 

p + m 



Sa+Ss \ Sb— Sa 
2m / m 



(3.10) 



(3.11) 



The functions ay and as are plotted in Figs. 2a and 2b for m = 1, 2, 4, 8 and cxd, the 
final curve being the bare propagator. From these plots we see that for large m, the 
deviation from the bare propagator due to the 1-loop self energy is dominated by the 
logarithmic contribution near the bare fermion mass shell, e = 0. With this in mind, 
we shall use the approximation 



Sb(x2 



Svr 4' 
Svr 4 



(3.12) 



(3.13) 



Taking e to be of order \/m for the purposes of the BS equation (see Eq. ( 3.16| ) below) 
these approximations give 



S{p) 



+ m 



m 



me 



— In - 



1 + 



In 771 



m 



(3.14) 



The vector and scalar parts of this approximate propagator (without the O (In m/m) 
corrections) are also plotted in Figs. 2a and 2b for comparison. 
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Turning now to the BS equation, we set the bound state momentum in Eq. ( |2.1| ) 
equal to = (2m + 6)iv^, where = (0,0, 1) and —6 is a "binding energy". This 
gives (according to the momentum distribution in Fig. 1) 



Tip)-- 
Setting 



(2vr) 



[m 



T{q)S 



[m + -)iv^ + 



in Eq. ( p. 14 ) gives 



S 



2 / 1. . Iql 

m \ 2 2m 



1 + 73 



+ O 



+ O 



Inm 



m 



7m- 
(3.15) 

(3.16) 



(3.17) 



Similarly 



S 



-{m + -)iv^ + qf. 



1-73 



(-l(5-ig3 + 



2m - 



— In — 



+ O 



Inm 



m 



(3.18) 



The |q| /2m term has been retained here to ensure convergence of the |q| integral in 
the BS equation below. Since the vertex T is defined with the fermion legs truncated, 
and 5 oc ^(1 ± 73), the only relevant part of T is the projection ^ (1 — 73) E ^ (1 + 73). 
With this in mind, the general forms in eqs. (|2.5D and (|2.6|) become. 



i(i-73)r^ni + 73)= Hi 



isj^gm, |q|) 



l(l-73)E^^l(l + 73)= ni-73)( 1 5(^73, |q|). 



(3.19) 



Substituting eqs. ( 3.17] ), ( p.l8| ) and ( 3.19 ) into Eq. ( p. 15 ) one obtains for the scalar 
states the single integral equation 



g{p) 



p.q 



9{q) 



(27r)3 {p-qf IpI |q| 



2 ' 



and for the axi-scalar states the single equation 

(fq 1 g{q) 



dip) 



(27r)3 (p-q)^ 



U + iqs+^-^-^ln-^ 



2m 4tt 2m v 2 



(3.20) 
2. (3.21) 



Note that, without the |q| /2m term in the denominator, translation invariance of the 
integrand implies that g is independent of |q|. In reality, g is a slowly varying function 
of |q|, and this extra 0(l/m^) term must be retained in e to account for the fact that 
the relevant region of integration in Eqs. ( |3.2C| ) and (3.21) extends out to 0{^/m) in 
the |q| direction, but only 0(1) in the q^ direction. Numerical solutions of Eqs. ( |3.20|) 



and (3.21) will be given in section 5. The function g is an even or odd function of q^ 
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corresponding to positronium states which are even or odd respectively under charge 
conjugation. 

In order to obtain a Schrodinger equation, we now rewrite the axi-scalar equation 
in the form 



9{P) 



(fq f"^ dq-i 1 

J-oo 2^ (p3-'73)2 + |p-q|' 



where 



(-i<5 + iq, + \^ + S+fe, |q|)) - iga + i + 5^-(93, |q 



.(3.22) 



(3.23) 



Assuming the integrand dies off sufficiently rapidly as gs — > —ioo, we deform the 
contour of integration around the pole at 



pole 
^3 = 



(3.24) 



to obtain 
dip) 



5(9^^^ iqi) 



(2^)' iP3 - + |P - qr -5 + + 25RS_(g3P°l^ |q|) ■ 



(3.25) 



(We could equally well deform the contour round the pole at (gf"'^)* if the integrand 
decays in the opposite direction, without affecting our final result.) 
Defining 

^(P3,|P|) = ^#^^^^^^^^ (3.26) 



<5 + ^ + 25RS_(p3,|p|) 



gives 



-5+^+2KS_(p3,|p|)U>(P3,|p| 
m 



Hqr'^'Acil). (3.27) 



(27r)2 (p3_^Pole)2 + |p_q|2 

In order to isolate the logarthmic infrared divergence we set 



pole , 
P3=P3 + 



(3.28) 



with fi small and real, and p^,"^'^ defined by analogy with Eq. ( p. 24] ). The right hand 
side of Eq. ( 3.27| ) then becomes 



r.h.s. 



(27r)2 ^2+o(^(|p|_|q|)) + |p_q|2mj 



F.T. of 



2n 



C + ln 



fir 

T 



(r) as /X ^ 0, 



(3.29) 
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where i?5>(q) 



|q|)- 



Following the reasoning of refs. |10| and [11|, this logarithmic divergence should 
be cancelled by the fermion self energy contribution 23ftS_(p3, |p|). However, from 
Eq. (3.23), we see that the logarithmic divergence in the self energy occurs at the bare 
fermion mass pole p^^'^'^ = —i(—5/2 + |p|^ /2rra), and not the dressed pole p^°^^- The 
problem lies in the use of the 1-loop approximation. If instead the fermion self energy 
is calculated to all orders in rainbow approximation, the self energy feeds back into 
the loop integral via the propagator to replace Eq. ( 3.23| ) by 



^-iP3, IpI) 




s-ip3 + - — h S-(P3, IpI 

2m 



(3.30) 



which provides a rainbow SD equation for $]_ in the non-relativistic limit. Then using 
Eqs. ( ^^) , (H), (l3^ ), (I3.29D and (|3.30D and fourier transforming we finally obtain 



|--V2 + ^(C7 + ln(mr))V 
y m 2tt J 



(3.31) 



agreeing with Eq. (|3.l| ). Had we started from the scalar equation ( |3.2[l| ) in place of the 
axi-scalar equation, the same result would have been obtained at Eq. ( p.29| ), leading 
to an identical Schrodinger equation. 

The important point to notice in this derivation is the significance of a non- 
perturbative solution to the SD equation in cancelling the infrared divergences. In 
the massless fermion limit, it is well known that chiral symmetry breaking plays a 
pivotal role in determining the bound state spectrum. It appears also that, even in 
the non-relativistic limit, the remnant effects of chiral symmetry breaking, via a non- 
perturbative solution to the SD equation, have a role to play. 



4 The Fermion Propagator 



The BSE described in Section 2 requires a fermion propagator input in the form of 
Eq. (|2.3| ) or Eq. ( |2.4D and this needs to be available over a region in the complex 
plane defined by Eq. ( p^ ) for 53 and |q| real. This is the region [13, |^ 



Q'^ =X + iY 



X > 



M2 



4 



(4.1) 



In this section we investigate ways of obtaining a solution to the SDE over 0. The 
fermion propagator, and thus the functions ay and as, must be well behaved over this 
region. 



The solution to the SDE (2^) is quite simple along the positive real (spacelike)-p^ 
axis. Substitution of the general expression for the fermion propagator ( |2.4D into the 
SDE gives an integral equation involving A and B functions which can be split into two 
coupled integral equations by simple projections. Angular integrals can be performed 
to leave one dimensional integrals (over the modulus of the q vector) , 



A(p2) - 1 



1 



47]-2p2 



dq 



qA{q^ 



p2 _|_ q2 



q^A^{q^) + B^{q^ 



4p 



In 



p + q 
p-q 
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The integrations range from to some UV cutoff along the positive real axis. This 
theory is super-renormalisable and thus has no ultraviolet divergences and so this 
cutoff is merely a numerical limit made large enough so that it has no bearing on the 
results. 

For a set of points p corresponding to the set of q points in the integration, the 
equations are iterated until convergence to leave the solution along the positive real 
axis. However, the solution is required for complex . We see three possibilities. The 
first is to use the converged functions A{q^) and B{q'^) in the integrals over the same 
contour (positive real q^) and supply the complex point p desired. The integrals should 
provide the solution at that point p. However the analytic structure of the integrands 
in Eq. ( [4.2| ) will not allow an analytic continuation by this method, because a pinch 
singularity in the integrand forces us to integrate through the point p 

The second possibility is to rotate the contour through an angle 2(/) in the p^ plane 
so that it passes through the desired point p |13]. In this way a cancellation of the 
complex parts within the logarithms occurs. Fig. 3 shows the first and second contours 
(Ci and C2 respectively). It can be seen from Eq. ( |4.2D that the logarithms will have 
real arguments along the radial portion of C2, while the arc portion contributes nothing 
to the integral because the integrand falls off sufficiently quickly in the ultraviolet |^. 

Based on the Landau gauge calculations of Maris |^, ^ we expect conjugate sin- 
gularities to occur in the second and third quadrants of the p^ plane away from the 
negative real (timelike) axis. Thus, as 2(j) increases towards tt from zero (and the neg- 
ative real p^ axis is approached) a singularity interferes and we may have convergence 
problems. It will be seen that these singularities can lie a fair way from (/> = ^ and 
convergence problems can occur for (j) not much more than ^ (ie: barely reaching into 
the second quadrant of p^). For the case m = to be considered shortly, no solution 
could be found for <j) greater than 0.90 radians (with a reasonable convergence criterion 
< 0.001). For this solution to be applied to the BSE we need to know the value of 
the fermion propagator for (j) from to ^ and so this method is not practical. However, 
although a slowing of convergence as (j) increases prevents a solution being attained in 
all of it does provide an accurate solution in a large portion of Q.. We therefore 
have a test for any A and B functions we wish to use in the BSE. 

The third possibility, and the one employed here and in ref. Q, is to find a good 
analytic fit along the positive real p^ axis and extend the solution into the complex 
plane by analytic continuation. These fits may be for A and B or for the functions 
ay and 175. The work of Maris suggests that it is not necessary for ay and 0$ to 
be entire functions for the fermions to be confined, only that there be no poles on the 
timelike p^ axis. Fits to functions A and B used in previous work Q based on the 
known asymptotic infrared and ultraviolet behaviour of these functions were tested by 
comparing them with the direct solution for various angles (j). The fits, adjusted to 
allow variable fermion mass m, are given by 

^fit(p^)= ^ /\ , +a3e-°^^^ + l, 
[ai + p^) 2 

Sfit(p') = + fese-'^P' + m. (4.3) 

02 +p^ 
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The parameters an,bn are functions of fermion mass. The numerical solution to which 
these functions were fitted is an iterative solution to the SDE using a non-uniform 51 
point grid along the positive real axis up to a momentum cutoff p = 1000 using a 0.1% 
tolerance in the integration routine. Plots of the numerical solutions and function fits 
for various m values are given in Fig. 4. 

Note that it is the a functions that are important in Eq. ( |2.1| ) and not A and 
B, and thus the effect of the fit on the denominator p'^A'^ + relating these must 
be considered. Conjugate poles exist where the factor p'^A'^ + B^ appearing in the 
denominator of the BSE integrand is zero. Table 1 lists the conjugate poles arising 
from the fits for each fermion mass and the corresponding maximum bound state 
masses allowed. The maximum M allowed is the value for which the boundary of Q in 
Eq. (4.1) coincides with the conjugate poles. No comment about the viability of our 
model BSE can be made until solutions are attempted because the integration region 
depends on the solution mass M. 

The location of the conjugate singularities for the m = case in Table 1 is slightly 
different to that reported in Ref. Q where it is —0.00400 it iO. 00666. This is because 
of the flexibility of the fitting functions. The fit in this work and that in Ref. |Q] for 
the zero fermion mass case had similar accuracy along the positive real p'^ axis but had 
the freedom to take on slightly different forms throughout the complex plane. This is 
because along the positive real axis the non-asymptotic form fixing parameters (04 
and 64) are only loosely determined. Despite the difference in the two results, the BSE 
calculation for bound state masses should show close agreement as each fit adequately 
models the direct solution throughout the complex plane. 

The singularities in the a fits for fermion masses greater than or equal to 0.1 lie on 
the negative real p^-axis. This suggests that free propagation occurs at these masses 
and the model is not confining. An accurate location of the singularities in the SDE 
solution would be needed before it can be said whether this result is due to the fits or 
the rainbow approximation used in the SDE solution. According to Ref. the rainbow 
approximation SDE solution is expected to be confining even for large fermion mass. 
Thus we assume our result is due to the lack of accuracy in our fits near the negative 
real p'^ axis and that it is likely that the singularities move close to that axis as m 
increases but never actually lie on that axis. 

Figs. 5a and 5b show plots of cry and as moduli respectively for zero fermion mass 
and angles (p = 0, (j) = ^ and (t) = \ against the p modulus (with a range far smaller 
than the UV cutoff used in our calculations). The direct solutions to the SDE and 
the fits are compared. It can be seen that the functions are very good fits along the 
positive real axis (0 = 0), where both ay and as are real. The fit is also good for 
(j) = ^. Real and imaginary components have not been given separately as they show 
similar agreement. In the case cj) = \ the fitting function has begun to deviate from 
the SDE solution. This is mostly due to the apparent difference in the location of a 
spike. Based on the largest bound state mass for m = reported in the next section, 
the BSE integration region extends along the direction = f out to a modulus 
of approximately 0.083. In this range the small angle solutions are very accurate but 
for larger (p, much of the error due to the difference in the location of the spike will 
be experienced. As the angle is increased further convergence problems occur until 
eventually no solution can be found at all [(j) > 0.90). 

The spike forming in these plots signals that, as (j) is increased, the contour of 
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integration approaches a singularity. In fact, the conjugate poles which lie just off (or 
on as is the case for larger m) the negative real axis (0 = ^) are approached. It is 
important that both the direct solution to the SDE and the fits used in this work have 
this feature. This spike was not seen in any other fits which we attempted. Based on 
Fig. 5 it seems clear that the direct solution to the SDE must have singularities close 
to those in the fitting functions. 

Because the spikes are not in exactly the same places some error will be introduced 
in the contributions from the large (j) part of $7. When the bound state mass becomes 
large, the large (j) contributions will become more important and thus we expect the 
error in the position of the spikes to result in some noise in the solutions to the BS 
equation for large fermion mass. 

The (T functions were studied for all fermion masses used in this work in the same 
fashion. The results were similar to the m = case and need not be shown here. In 
each case, when (j) was increased far enough, a spike was observed in both the fit and 
solution, after which lack of convergence prevented an SDE solution. 

However, for very large fermion masses, the accuracy of the fits decreases as m 
increases, and with good reason. As m tends to infinity, the functions A and B 
approach constants (1 and m respectively). For moderately large fermion masses 
experienced in this work, these functions become almost constant along the positive 
real axis while having a singularity near the negative real p^-axis. It is too much to 
ask for simple four parameter fits along the positive real p^-axis to reproduce accurately 
complex behaviour deep into the real timelike p^-axis. The 1-loop propagators. Fig. 2 
described in section 3 illustrate this well. There one can see how smooth and level the 
a functions are along the positive real axis and also how steep the functions become 
back along the negative real axis. 

Before moving on to the next section, we return briefly to the 1-loop approximation 
to the fermion propagator necessary for the non-relativistic approximations described 
in section 3. Fig. 6a compares our rainbow approximation solution A to the 1-loop 
result given in Eqs. (^^) and (|3.4| ). Fig. 6b compares B from our rainbow approxima- 
tion solution and the result in Eqs. ( f3.5| ) and ( |3.6[ ). Both of these comparisons were 
made at a large fermion mass {m = 5). It can be seen that the curves in each case are 
in reasonable agreement, at least for spacelike momenta. 



5 Numerical Solution of the Bethe-Salpeter Equation 



The fits given by Eq. ( |4.3| ) to the fermion propagator for a range of fermion masses 
were used in the solution of the Bethe-Salpeter coupled integral equations Eq. (|2.7|) . 
This problem was restated in Eq. as an eigenvalue problem. 

A grid of 25 x 25 (|q|,Q'3) tiles were used for the iterative procedure with linear 
interpolation on each of those tiles used for the sums {Tijfj) which are supplied at 
the corners of the tiles from the previous iteration. The tiles were non-uniform in 
size and an upper limit to the momentum components (|q| and q^) of between 3.0 
and 9.0 was used. The equations were iterated to convergence each time to determine 
eigenvalues for a given test bound state mass M. The Bound state masses were located 
by repetitive linear interpolation or extrapolation to search for the point where the 
eigenvalue A of Eq. ( 2.131) is 1. This was repeated for each of the fermion masses 



ranging from to 5.0. This procedure was used for each of the four non-degenerate 
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bound state symmetries described in the appendix. 

Table 2 shows the bound state masses for each of the four symmetries (scalar 
C = +1, scalar C = — 1, axi-scalar C = +1 and axi-scalar C = — 1) for all fermion 
masses considered. Fig. 7a displays the solutions M for fermion mass 0-0.1. Fig. 7b 
shows M — 2m over the greater range of 0-5. The axi-scalar C = +1 solution is a 
degenerate axi-scalar /axi-pseudoscalar pair of Goldstone bosons for the case m = 0, 
as seen in previous work Q. Minor differences between Ref. [Q| and the current work 
at m = are due to small differences in the propagator fits, as explained in section 4. 

For small m the bound state masses rise rapidly with with increasing fermion mass. 
The mass of the "Goldstone" axi-scalar C = +1 state scales roughly with the square 
root of the fermion mass, in agreement with the Gell-Mann-Okubo mass formula 
In fact, for fermion masses to 0.1 a linear regression against y/m has correlation 
coefficient 0.9964 with the mass growing as approximately 1.27 x -^/m. (The accuracy of 
the solution at m = 0.001, which comes out with an anomalously low bound state mass, 
is severely affected by numerical inaccuracy arising from the sensitivity the bound state 
mass to the eigenvalue A in Eq. ( p.l3| ).) For large fermion masses, the bound state mass 
rises predominantly as twice the fermion mass plus possible logarithmic corrections. 
However, there appears to be a good deal of noise in the large m solutions, reflecting 
the difficulty in accurately modelling the fermion propagator deep into the timelike 
region from spacelike fits. No solutions corresponding to states of negative charge 
parity were found for m > 1.0. 

Numerical solutions to the integral equations ( |3.20| ), ( |3.21 ) arising from our non- 
relativistic treatment are listed in Table 3 and plotted in Fig. 7b. Solutions with 
positive 6 were found for fermion masses m > 1.0 in the positive charge parity sector. 
We were unable to locate any solutions to Eqs. ( [3.20 ), ( 3.21 ) corresponding to negative 
charge parity states over a broad range of 5. Also given in Table 3 and Fig. 7b are 
the two lowest lying s-wave solutions to the Schrodinger equation from the numerical 
work of Tarn et al. [^, given by Eq. ( |3.2D . 

The lack of exact agreement between the non-relativistic, 1-loop approximations 
Eqs. ( p. 20 ) and ( 3.21| ), and the Schrodinger equation result Eq. (3^) is to be expected. 
As pointed out in Section 3, a complete cancellation of infrared divergences can only 
occur if the fermion self energy is calculated non-perturbatively to all orders. From 
Table 3, we see that at very high fermion masses, the accuracy of the 1-loop approxi- 
mation is significantly affected as the conjugate poles in the propagator, measured in 
momenta scaled by the fermion mass, move closer to the bare fermion mass pole (see 
Fig. 2). At more moderate fermion masses, m ~ 5, the 1-loop approximation is more 
respectable. 

We see no clear agreement between the numerical results of Eq. (2^), and either 



non-relativistic approximation Eqs. (|3.20| ), (3.21), or the Schrodinger equation result 
Eq. ( |3.2D . Our analysis of the non-relativistic limit of the BS equation exposes the 
importance of the analytic structure of the fermion propagator in the vicinity of the 
bare fermion mass pole p'^ = —m?. The uneven nature of the lower two curves in Fig. 7b 
indicates that the determination of the timelike fermion propagator by an analytic fit 
to the spacelike propagator is inadequate for fermion masses m > 1. It is clear that 
a more careful analysis of the timelike nature of the fermion propagator, possibly 
involving a fully non-perturbative treatment of the SD equation to include remnant 
chiral symmetry breaking, is necessary for determining the bound state spectrum for 
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even moderately large fermion masses. 

It is important to note that the poles in the fermion propagator fits listed in Table 1 
lie outside the BS integration region $7 for all solutions obtained. This can be verified 
by observing that all masses in Table 2 are lower than the values Mmax listed in Table 1. 
A similar situation arises for the non-relativistic limit calculations. Listed in Table 3 
are maximum allowed 6 values if the integration region sampled by Eqs. ( |3.2[1| ) and 
( 3.2l|) is not to impinge on the conjugate propagator poles ^3°^'^ and (^3°^'^)* defined 
in Eq. ( |3.24| ). In all cases the numerical results lie within the permitted region. This 
requirement is equivalent to demanding that §3°'^ should not cross the real axis as 
as |q| ranges from to cxd. Interestingly, such a crossing would entail a more careful 
evaluation of residues than that carried out in Section 3 leading to the Schrodinger 
equation. 

We note that the Schrodinger equation results of Ref. include the first five 
s-wave states. It would certainly be of interest to locate the excited states within 
the framework of our BS treatment of QED3. We have searched for solutions to 
the eigenvalue equation ( p. 13 ) corresponding to excited states, and find in general 
no solutions within the mass ranges allowed by the values M^ax in Table 1. Since 
there is no reason to assume that the s-wave spectrum should be bounded above, it 
seems likely that there will be solutions to the BS equation for which the region of 
integration Q does include the conjugate propagator poles discussed in Section 4. It 
follows that the functions /, C/, V and W in the BS amplitudes of these states should 
have compensating zeros, in order that the right hand side of the BS equation be 
integrable. We conjecture that, if the fermion propagator has an infinite set of poles, 
there will be a sequence of excited states, the nth excited state having n pairs of zeros 
in its BS amplitude. This conjecture is consistent with the the first excited state of the 
Schrodinger equation, also listed in Table 3, for which the wave function has a single 
zero. 

Although we are unable to determine accurately the spectrum in the large fermion 
mass limit, our calculations strongly suggest that there are no scalar or axi-scalar states 
with negative charge parity in this limit. This is consistent with the non-relativistic 
quark model in four dimensions in which negative charge parity scalar and pseudoscalar 
states are forbidden by the generalised Pauli exclusion principle [14|. We note, however, 
that there is nothing to exclude such states in a fully relativistic BS treatment 
and indeed, negative charge parity scalar and axi-scalar states are found within the 
current model for light fermions. 



6 Conclusions 

In this paper we have solved the combination of rainbow Schwinger-Dyson and ho- 
mogeneous Bethe-Salpeter equations in the quenched ladder approximation for three 
dimensional QED with massive fermions. QED3 was chosen because, like QCD, it is 
confining but without the complications of being non-abelian. A four-component ver- 
sion of this theory is used because, also like QCD, it provides a parity invariant action 
with a spontaneously broken chiral-like symmetry in the massless limit. The approxi- 
mation is amenable to numerical solution, and should help assess the limitations of a 
technique frequently employed in models of QCD [^. 

The work in this paper carries on from a previous study of the same subject but 
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with the fohowing extensions. Firstly, non-zero fermion masses is considered. Secondly, 
an analysis of the fermion propagator in the complex plane is carried out in order to 
assess the appropriateness of the approximations involved. Thirdly, an analysis of the 
non-relativistic limit, i.e., large bare fermion mass, is made in an attempt to compare 
with existing Schrodinger equation studies of QED3. 

The rainbow SD equation was solved in Euclidean space to give a fermion propaga- 
tor for spacelike momenta, Euclidean > 0. The propagator is chirally asymmetric, 
and in the massless fermion limit, gives rise to a doublet of massless Goldstone positro- 
nium states analogous to the pion. Solution of the BS equation for massive positronium 
states requires knowledge of the fermion propagator S (p) in the complex p^-plane ex- 
tending away from the spacelike axis, and a finite distance into the timelike axis < 0. 
By rotating the contour of integration we were able to extend the spacelike solution 
into part of the complex plane. However, the occurrence of complex conjugate poles in 
the fermion propagator prevented a numerical solution to the SD equation throughout 
the complete region of the complex plane sampled by the BS equation. This forced us 
to apply analytic fits to the propagator along the positive real p^-axis for use over the 
required part of the complex plane. 

Our propagator fits were found to have conjugate poles located close to those of the 
direct solution for small to moderate fermion masses. This, combined with the accuracy 
of the fits throughout much the complex plane, made our choice of propagator very 
attractive. The singularities in the fits were found to move onto the negative real p^- 
the fermion mass increased. This was not interpreted as a loss of confinement 
but instead attributed to a lack of accuracy in the fits deep into the timelike region as 
the fermion mass became large. This reduction in accuracy of the fits for large m was 
due to the nature of the functions along the positive real p^-axis where the fits were 
made, and the presence of a singularity near the negative real p-^-axis in the vicinity 
of the bare fermion mass pole p'^ = —m?, but off the timelike axis. 

BS solutions were found for four pairs of parity degenerate states. These pairs 
were the scalar /pseudoscalar C = +1 and C = — 1 and the axi-scalar/axi-pseudoscalar 
C = +1 and C = — 1 states. For small to moderate fermion mass the bound state mass 
was found to increase smoothly with m. The axi-scalar C=+l doublet, analogous to 
the pion, was the lowest in energy, with a mass rising roughly with the square root 
of the bare fermion mass. For moderately large bare fermion masses {m/e^ greater 
than unity) the positronium masses rise as twice the bare fermion mass, plus a possible 
logarithmic correction. However, an unacceptable level of noise was found to develop in 
our results for these larger masses, which we attribute to inaccuracies in the analytically 
continued fermion propagators in the important region near the bare fermion mass pole. 
No negative charge parity (C = —1) solutions were found for bare fermion masses above 
m/e^ ~ 1.0, consistent with the generalised Pauli exclusion principle of non-relativistic 
QCD4. 

The conjugate poles in the fermion propagators were found to keep clear of the 
integration regions required for the BS solutions for the lowest state in each of the 
four space parity/charge parity sectors considered. However, it appeared that this 
would not be so for any excited states. We therefore conjecture that the excited 
positronium states have zeros in their BS amplitudes positioned so as to cancel the 
poles in the propagators encountered within the integral in the BS equation ( |2.1| ). This 
requirement of compensating zeros was too demanding on our current numerical code. 
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and as a result, no excited states were found. 

In vector calculations under way at present, where the bound state masses are 
expected to be larger, the conjugate poles in the fermion propagator seen in this 
work may interfere. Since the fits used in this work appear to have their singularities 
close to those in the actual Schwinger-Dyson solution, we may find that the rainbow 
approximation and the resulting propagator fits will be inadequate for a study of 
vector states in QED3. This is a very challenging problem and we hope to report on 
our results in the near future. 

A non-relativistic analysis of the BS equation was also carried out assuming, in 
the first instance, a 1-loop approximation to the fermion propagator. However it was 
shown that, in order to cancel infrared divergences completely between the photon 
propagator and fermion self energy, as proposed by Sen ||l^ and Cornwall it is 
necessary to evaluate the fermion self energy non-perturbatively. Only if this is done 
can the Schrodinger equation be rigorously obtained in the large fermion mass limit. 
In spite of this, numerical solutions of the 1-loop equations give reasonable agreement 
with the Schrodinger equation for moderately large fermion masses m/e^ ~ 5. 

In summary, we were able to carry out an acceptable analysis of the bound state 
spectrum of QED3 near the chiral limit m — > by using analytic fits to the spacelike 
fermion propagators in the BS Bethe-Salpeter equation, and in the non-relativistic limit 
m — > oo by expanding to lowest order in inverse powers of the fermion mass to obtain 
a Schrodinger equation. However, there remains an intermediate mass range m/e^ ~ 1 
for which neither of these techniques is adequate. It is clear that a more careful non- 
perturbative analysis of the fermion propagator in the vicinity of the bare fermion mass 
pole is necessary before an accurate determination of the QED3 positronium spectrum 
at intermediate fermion masses can be made. If a direct analogy with QCD models 
based on the Bethe-Salpeter equations is made, we conclude that particular care must 
be taken in modelling quark propagators for quarks whose mass is close to the mass 
scale of the theory, namely charm quarks. 

Appendix - Transformation Properties in QED3 

The four-component QED3 action in Minkowski space |jl6| 



involves 4x4 matrices 7^ which satisfy {7/^,7;/} = 2?7^j,y where rj^u = diag(l, —1, —1) 
with ^ = 0, 1 and 2. These three matrices belong to a complete set of 16 matrices 
{ia} = 74, 75, 745, 7m, 7^4:7/^5, 7/^45} satisfying ltr(7A7'^) = ^a'^ 



71,2 



70 = 


[ 0-3 






C 


;) 


7m5 






4 / u J 1 5/0 \ 45 

74 = 7 = I r n I ' '^^^'^ = \ il J ' '^'^^^'^ " -^475, 

u4,/i5or/i45 _ „/i!^^ ^ ^ 



The three 7^, and 74 and 75 are five mutually anti-commuting matrices. This is 
unlike the 4-dimensional case where no analogue of 74 exists. 
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The action Eq. ( [A.l[ ) in the massless case m = exhibits global U{2) symmetry with 
generators {1,74,75,745} which is broken by the generation of a dynamical fermion 
mass ||l7| , |l^ to a C/(l) x U{1) symmetry {1,745}. The action is also invariant with 
respect to discrete parity and charge conjugation symmetries, which for the fermion 
fields are given by 

ip{x) ip'{x') = nV'Cx), ^{x) i^\x') = v^(x)n-\ (A.2) 

V'(x) ^ ij'{x) = C^ix)^, i^{x) i^\x) = -i){x)^C\ (A.3) 

where x' = (x", — x^, x^). The matrices 11 and C are each determined only up to an 
arbitrary phase by the condition that the action Eq. ( [A.l| ) be invariant 

n = 7i4e*'^^^45 ^ c = 72e*'^^^« , (0 < (/>p, 0c < 27r) (A.4) 

Scalars, pseudoscalars, axi-scalars and axi-pseudoscalars are defined by the follow- 
ing transformation properties under parity transformations 



$^(x) - 


^^'{x') = 


^P'{x) - 


^P^'ix')-- 








^^P^'{x') 



= i?P$^^(x), 

= -Rp^^P^{x), (A.5) 
where <i>^'^ and (^^^^ are doublet states <I> = (<I'4, <1*5)^, and 

\ —sm Z0P cosz<pp / 

Similar transformation properties exist for charge conjugation. 

The most general forms of the Bethe-Salpeter amplitudes |15| for bound scalar and 
pseudoscalar states are 

T^{q,P) = If + iig+ n + e^upP^'q^l'^'k, (A.7) 
r^^(<?,P) = ^45r'{q,P), (A.8) 

where /, g, h and k are functions only of q^, and q ■ P. BS amplitudes corresponding 
to the components ^>4 and <I>5 of axi-scalars and axi-pseudoscalars take the general 
form 

and 

[r^'\,,P)) =^''[r^'H,,P)) ■ ^^-'"^ 

Furthermore, the charge parity C = ±1 of the bound states is determined by the 
parity of the functions /, g, h and k under the transformation q ■ P ^ —q ■ P. The 
quantity q-P is the only Lorentz invariant which changes sign under charge conjugation 
and thus determines the charge parity of those functions. 

Our conventions for Euclidean space quantities are summarised in Appendix A of 
Ref. [Q]. In particular Euclidean momenta and Dirac matrices are defined by 

p(E) _ _ . p(M) p(E) _ p(M) (E) _ (M) (E) _ (M) 
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m 






0.000 


-0.0034 ±i 0.0057 


0.142 


0.001 


-0.0041 ±i 0.0064 


0.153 


0.004 


-0.0060 ±i 0.0086 


0.182 


0.009 


-0.0081 ±i 0.0140 


0.206 


0.016 


-0.0121 ±i 0.0192 


0.247 


0.025 


-0.0216 ±/ 0.02(30 


0.325 


0.036 


-0.0314 ±i 0.0345 


0.386 


0.049 


-0.0468 ±i 0.0387 


0.464 


0.064 


-0.0618 ±i 0.0417 


0.522 


0.081 


-0.0815 ±i 0.0440 


0.590 


0.1 


-0.0647 ±i 0.0000 


0.509 


0.5 


-0.4894 ±i 0.0000 


1.399 


1 


-1.4260 ±i 0.0000 


2.388 


2 


-4.8925 ±i 0.0000 


4.424 


3 


-10.3776 ±i 0.0000 


6.443 


4 


-17.9613 ±i 0.0000 


8.476 


5 


-27.3616 ±i 0.0000 


10.462 



Table 1: Conjugate singularities for fermion propagator fit and corresponding limits 
on bound state mass for fermion masses from to 5.0. 
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m 


Scalar C = +1 


Scalar C = — 1 


Axi-scalar C = +1 


Axi-scalar C = — 1 


[Ref. i] 


0.080 ± 0.001 


0.123 ± 0.002 





0.111 ± 0.002 





0.077 


0.118 





0.108 


0.001 


0.087 


0.126 


0.004 


0.116 


0.004 


0.110 


0.151 


0.054 


0.140 


0.009 


0.140 


0.178 


0.090 


0.167 


0.016 


0.175 


0.217 


0.127 


0.204 


0.025 


0.215 


0.269 


0.167 


0.254 


0.036 


0.256 


0.316 


0.208 


0.300 


0.049 


0.298 


0.367 


0.248 


0.350 


0.064 


0.343 


0.411 


0.293 


0.390 


0.081 


0.389 


0.456 


0.340 


0.431 


0.1 


0.439 


0.496 


0.391 


0.479 


0.5 


1.311 


1.388 


1.261 


1.352 


1 


2.297 


2.387 


2.243 


2.336 


2 


4.330 




4.233 




3 


6.348 




6.227 




4 


8.379 




8.243 




5 


10.365 




10.219 





Table 2: Bound state masses for fermion masses from to 5.0. (All masses ±0.001 
unless otherwise stated). The Axi-scalar C = +1 solution with m = stated here is 
an analytic result. 
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m 


^max 


Scalar 
Eq.(^.20|) 


Axi-scalar 
Eq.(|3.21) 


Eq.(^ with 
A = Ao 


Eq.(|33) with 
A = Ai 


1 


0.332 


0.285 


0.262 


0.322 


0.503 


2 


0.421 


0.371 


0.338 


0.377 


0.558 


3 


0.473 


0.419 


0.381 


0.409 


0.590 


4 


0.511 


0.452 


0.410 


0.432 


0.613 


5 


0.540 


0.476 


0.433 


0.450 


0.631 


100 


0.947 


0.792 


0.734 


0.688 


0.869 


1000 


1.272 


1.030 


0.968 


0.872 


1.052 



Table 3: Non-relativistic 5 solutions for positive charge parity, from Eqs. (3.20), 
( 3.211) and the two lightest s-wave Schrodinger equation results of Tarn et al. 
Eq. ( p.2| ) using Aq = 1.7969 and Ai = 2.9316. The first column contains the maxi- 
mum values of 6 allowed before conjugate singularities arise in the fermion propa- 
gator used in our non-relativistic calculations. 



23 



Figures 



Figure 1: Diagrammatic representation of Eq. (p7 



Figure 2: Figures 2a and 2b show 1-loop approximations using eqs.(2^), {2A) and 
( |3.3D - ( |3.6D for rn^av and mas respectively (solid lines). These are compared 
with the vector and scalar parts of the approximation Eq. ( |3.14D (dashed lines) . 
The curves are drawn for fermion masses (from bottom to top) 1, 2, 4, 8, and oo. 



Figure 3: The first and second (deformed) contours of integration Ci and C2 for 
solution to Eq. (|4^. 

Figure 4: Figure 4a compares the SDE solutions and fitting functions for A — 1 for 
fermion masses (from top to bottom) m =0, 0.025, 0.1, 1 and 5. Figure 4b shows 
B — m for fermion masses (from bottom to top) m =0, 0.025, 0.1, 1 and 5. 

Figure 5: Figures 5a and 5b show SDE solutions and function fits for ay and as 
respectively for fermion mass with angles = (O), ^ (+) and ^ (□). 

Figure 6: Figure 6a compares the function Aij)^) from the rainbow SDE calculation 
(solid curve) and the 1-loop result (dashed curve) and figure 6b compares -B(p^) 
results along the positive real p axis for fermion mass m = 5.0. 

Figure 7: Figure 7a shows bound state masses (M) against fermion mass 0-0.1. Fig- 
ure 7b is a plot of M — 2m for m =0-5. In each plot the scalar C = +1 (O), scalar 
C = —1 (-|-), axi-scalar C = +1 (□) and axi-scalar C = —1 (x) states are drawn 
with solid curves. The non-relativistic predictions of Eq. ( |3.20D and Eq. ( |3.21| ) 
are the scalar C = +1 (O) and axi-scalar C = +1 (□) states respectively and 
are drawn with dashed lines. Eq. ( p.2[ ) with A = Aq (lower solid curve with no 
symbols) and with Ai (upper solid curve with no symbols) are also plotted in 
figure 7b. 
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1.011 
1.01 
1.009 
1.008 
1.007 
1.006 
1.005 
1.004 




Fermion mass (m) 



